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Abstract 

The first excitation energy in the two-level pairing model is investigated in terms 
of the equilibrium and the small fluctuation around it. The basic idea is an extension of 
results presented in a previous paper by the present authors. In this investigation, the 
result obtained in the previous paper plays a central role: At the limit of the weak and 
the strong interaction strength, the results are in good agreement with the exact one. 
The former and the latter are calculated in the framework of the su(2) <S> su(l, 1)- and 
the su(2) <S> sw(2)-coherent state, respectively. On the basis of the above conclusion, 
the intermediate region for the interaction strength is described in terms of the idea of 
the interpolation and it is shown that the agreement of the result with the exact one 
is quite good. 
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The description of the first excited state in the two-level pairing model for the closed shell 
system has been regarded as a basic, but a modest problem in many-body physics. For exam- 
ple, the recent investigations can be found in Ref. 1). Following this quiet stream, recently, 
the present authors have considered this problem from tow slightly different viewpoints. One 
was reported in Ref. 2) and the other in Refs.3) and 4). Especially, in Ref. 4) (referred to 
as (I)), the ground-state and the first excitation energy are calculated in the framework of 
three different coherent states which consisit of four kinds of boson operators. In (I), they are 
called (i) the su(2)®su(l, 1)-, (ii) the su(l, 1)- and (iii) the su(2) <S> su(2)-coherent 

states. The picture adopted in (I) is based on the equilibrium and the fluctuation around 
it and the equilibrium state is described by a coherent state. Therefore, we obtain three 
different equilibrium states. Roughly speaking, for the first excited energy, in the region of 
weak strength of the pairing interaction, the state (i) gives rather good result. But, in the 
region of strong strength, the state (iii) gives good result. Compared with the above two 
states, it seems to us that the state (ii) gives a slightly worse result. Of course, the above 
rough summary is based on the comparison with the exact result. 

Before entering the central parts of this paper, we mention our basic viewpoint. We pay 
attention to the relation appearing in the appendix (B) of (I), i.e., (I-B-l). This relation 
demonstrates that, if we solve our system exactly in the quantum mechanical framework, 
the result does not depend on the choice of the function g(K). The part of the Hamiltonian, 
Hi, is recast into 



H x = -G 
= -G 



g(K) ■ Vhc*\/2T + K\J 2M - K\J 2L - K ■ (^J g(K + h)^ + h.c. 



Vhc*\ 2T + K\ 2M - K\ 2L-K + h.c. 



where h.c. means the Hermitian conjugation. It should be noted that the relation holds 
for any form of g{K) if it can be defined. On the other hand, the c- number replacement 
gives the form 



Hi = -G 
= -G 



y/g(K) ■ VhcW2T + KV2M - K^2L - K ■ [y/g(K + K)J + c.c. 
\/ftcV2T + KV2M - KV2L - K ■ y/g(K) ■ g(K + h)^ + c.c.l , (2) 



where c.c. means the complex conjugation. The relation (T5]) tells us that the classical 
solution based on the relation (J2j) depends on the choice of g{K). The above means that the 
classical counterpart under the c-number replacement cannot be fixed uniquely, that is, there 
exist infinite possibilities for the classical solutions which should correspond to the unique 
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quantum mechanical solution. However, we must note that we determine the equilibrium by 
a classical solution and the fluctuation is given in the lowest order for classical and quantum 
mechanical framework. Therefore, the result depends on the form of g(K), and the principle, 
with the aid of which g(K) is fixed, is required. In (I), we investigated three possibilities 
which come from three forms of the coherent state. 

As was demonstrated in (I), the key to obtain a good result is to choose the function 
f(K) related to g(K) and defined through the form (I-B-2), i.e., f(K). The three cases are 
shown in the relation (1-4-4) and in Fig. (1-2), the different behaviors are explicitly presented. 
Here, "good result" means the agreement with the exact one. From this figure, we learn 
that if we succeed in finding the function f(K), which satisfies the following condition, the 
good result may be expected: For a small value of K, f(K) shows the behavior shown in 
Fig.(I-2a) and in the region K^Q/2, the behavior of f(K) should be as shown in Fig.(I-2c). 
In this paper, we present some candidates for f(K) which show the above behavior under 
the idea of an interpolation. 

Instead of f{K), we investigate the function n(K) defined as 

f(K) = 2G»(K) . (3) 

In the case of the states (i) and (iii), fi(K) can be expressed in the form 

m(K) = y/K(K + l)(Q - K) , (4) 
(i m (K) = K(Q - K) . (5) 

Hereafter, we use the unit K — 1. The function fii(K) can be expanded in the form 

m(K) = Vk(Q + (Q/2-1)K + ---) . (6) 

On the other hand, fim(K) can be rewritten as 

fi m (K) = (f2/2) 2 + x {Q/2 - K) - {(2/2 - Kf . (7) 

Therefore, the leading terms for miK) near K ~ and for n\w(K) near K ~ Q/2 are as 
follows: 

m(K) = Vkq , (8) 

fi m {K) = (Q/2) 2 + x (Q/2 - K) . (9) 

Our problem is to find an explicit form for fi(K) which reduces to the forms jS]) and (Q for 
K ~ and K ~ Q/2, respectively, in the framework of the idea of an interpolation. 
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First, let us search fi(K) in the form 

H(K) = s/K(K + n)A(B - K) . 



(10) 



Here, n, A and B denote parameters to be determined as functions of Q. The function fi(K) 
is expanded in the form 



'Ha) 



fi{K) = VK[VnAB + A{B/2^- y/n)K + • • •] . 
Also, in terms of (f2/2 — K), fi(K) is expressed as 



fx(K) = ^Q{Q + 2n)/2 ■ A{B - Q/2) 



+ ^Q{Q + 2n)/2 ■ A 
A 



i - (B-n/2) 



2(Q + n) 
J?(J? + 2n) 



V / i7(i7 + 2n) _ 



{p + n) + (B- Q/2) 



n 



Q(Q + 2n)_ 



(Q/2-K) 
{Q/2 -Kf + - 



lib) 



Comparison of the forms filial) and (lllbl) with the asymptotic forms (JE|) and (J9]) gives us the 
relation 



^AB = Q , 

^Q{Q + 2n)A(B - Q/2) = (Q/2f 
2(Q + n) 



1-(B- Q/2) 



Q(Q + 2n) 



(12a) 
(12b) 

(12c) 



The relations ( Il2bl) and ( 11 2cf) lead to 

_ VQ(Q + n) 



,3 ' 



B = Q 



Q + 3n/2\ 



(y/n + 2n) 

Substituting the form f|T3l) into the relation f)12al) . we have 

(Q + 2nf 



Q + n 



) 



(13) 



n 



Q(Q + 3n/2) 2 



By solving the relation (|14p . we can express n in terms of Q and the relation (fT5|) ^ 
and B expressed in terms of Q. In the case Q = 5, for which we show the numerical 
in figures, we have 



(14) 

;ives A 
results 



n = 2.8441 



A = 0.5020 



B = 5.9064 



(15) 
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Next, we investigate a possible improvement of the form (FlOi) with the relations ( fl3i) and 
dHj). We set up the following form: 



p(K) = \JK{K + n)A(B - K) ■ p(K) 



(16) 



If p(K) = 1 in the region < K < ft/2, the improvement is not achieved. If we intend to 
keep the forms f)13p and (I14p in the improvement, p(K) should satisfy 



p(0) = p(ft/2) = 1 , p'(ft/2) = 
As a possible candidate, we can choose the form 



(17) 



p(K) 



l + C ■ (ft/2)-\ft/2- Kf 
l + C ■ {ft/2)- 2 {ft/2- K) 2 ' 



(18) 



Here, C denotes a new parameter to be determined as a function of ft. From the comparison 
with the exact result, we can see that the large value of K (~ ft/2) should be improved. 
Then, we expand p(K) in terms of (ft/2 — K)\ 



p(K) = 1 - (C/2) • (ft/2)- 2 (ft/2 — K) 2 + ■ 
Then, the expression (lllbl) should be improved in the form 



(19) 



the term improved from the third term in the expansion flllbj) 



.4 



v / ft(ft + 2n) 



(ft + n) + (B-ft/2) 



n 



ft{ft + 2n) 



+ y/ft(ft + 2n)/2 ■ A(B - ft/2) ■ (C/2)(ft/2) 



-2 



[ft/2 - Kf . 



(20) 



By putting the form ( 1201) equal to the third terms of the relation (J7j), we obtain 



.4 



^ft(ft + 2n) 



(ft + n) + (B- (2/2)- 



n 



ft(ft + 2n) 

+ ^ft(ft + 2n)/2 ■ A(B - ft/2) ■ (C/2){ft/2)- 2 ■- 



1 . 



By substituting the result (IT51) for A and B into the form (|2~TI) . C is determined as 

2n(ft + 5n/4) 



C 



(ft + 2n) 2 



(21) 



(22) 



The case (2 = 5 gives us 



C = 0.4260 



(23) 
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As is clear from the relation (I-B-4), it is also an important task for the present description 
to obtain the function g(K) defined through the relation 



H(K) = y/K(K + l)(n - K)^g(K)-g(K+l)- 1 
In the present case, we have 

For the expression ffTUj) . we have 

Zi{n,K) = i. 

For the expression ffTB"]) with the candidate ([TBI . K) is given as 

£ 2 (f2, K) = {Q/2)- 2 \{Q/2) {(1 + 1/V2\/C) + i ■ 1/V2 <fc) - K 



(24) 



(25) 



(26) 



{Q/2) {(1 - l/V2</C) + i ■ l/V2<fc) - K 



x 



(n/2)(i + i/vc)-K 



-1 



x [c.c] . 



(27) 



The notation [c.c] denotes the complex conjugate of all the previous terms. Hereafter, we 
distinguish the two cases by the indeces 1 and 2. For finding the functions gi{K) and g2(K), 
the following formulae are useful: 



(i) if h(K) -h(K + iy l 

(ii) if h{K) ■ h^K+iy 1 
(m) if h(K) ■ h(K + I)- 1 



a , h(K) = (3a~ K , 

a - K , h(K) = pr(a + 1-K) 

a + K, h(K) = (3r(a + K)- 1 . 



(28a) 
(28b) 
(28c) 



Here, a and (3 do not depend on K and r(z) denotes the gamma-function. As can be seen 
in the relation (I-B-4), we need only the form g'(Ko)/g(Ko), and then, it may be enough to 
put (3=1 for the formula (|28|) . With the use of the formula (1281) . we obtain the following 
relations for g{K): 

1 2 



g x {K) = A~ 2K ■ r(l + K) ■ r(n + K)~ l ■ [r(B + 1 - K) ■ r(Q + 1 - K)- 1 ] ' 
g 2 {K)= 9l {K){Q/2f K 

xr((/2/2){(l + \/^2<fC) + i • l/V2<fC} + 1-K) 

xr((n/2){(i - i/V2~Vc) + % ■ i/V2~Vc} + i-k) 



(29) 



x 



r((f2/2)(l+i/VC) + l- K) 



-1 



x [c.c] 



(30) 
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Fig. 1. The energies obtained by the use of the various states are depicted as a function of G. 
Here, En\ and ifym denote the derived energies by using of m and in Eqs.((l]) and ©, 
respectively, and E\ denotes the energy by using of \i{K) in (fT0j) with ([lljjW ([T5j) . where Q = 5. 
Here, E exSuCt shows the exact eigenvalue of the Hamiltonian. 
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G 

Fig. 2. The frequencies are depicted as a function of G. Here, un\ and uau) denote the derived 
frequencies by using of [i\ and fj,m in Eqs.(|4|) and (|5|), respectively, and u\ denotes the frequency 
by using of fi(K) in (fTUj) with (fT3lW(fT5j) . where ft = 5. Here, u ex shows the exact result. 

In FigJU the derived energies of equilibrium are compared with the exact ground state 
energy E exact . Here, and denote the energies obtained by using /z; in Eq.Q and 
yUjii in Eq.([5]), respectively, and E\ denotes the energy newly obtained by using the \i{K) in 
Eq. fflO]) . Also, in FigfSl the derived first excitation energies or frequencies of small oscillations 
around the equilibrium are compared with the exact first excitation energy u ex . Here, uim 
and o;(iii) denote the first excitation energies obtained by using /i; in Eq.(jl]) and in Eq.(JS]), 
respectively, and lo\ denotes the first excitation energy newly obtained by using the fi(K) in 
Eq. (fTUj) . The energy of equilibrium is not so refined except for the region presenting a dip 
structure in -E(iii). However, the frequency Ui is obviously refined. In the region of the weak 
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Fig. 3. The energies obtained by the use of the new /x(-PC) are depicted as a function of G. Here, 
Ei and E^ denote the derived energies by using of p(K) in (fT0|) and (fl"6j) with (fT8|). respectively, 
and inexact shows the exact eigenvalue of the Hamiltonian. Here, Q = 5. 
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Fig. 4. The frequencies obtained by the use of the new p(K) are depicted as a function of G. 
Here, w\ and W2 denote the derived frequencies by using of p>(K) in (fl~Qj) and (fl~6j) with (fTHl) . 
respectively, and lu cx shows the exact result. Here, 17 = 5. 

interaction strength G, the behavior of oj\ is similar to uj^ by using the su(2) <g) su(l, 1)- 
coherent state as is expected. On the other hand, in the region of the strong interaction 
strength, the behavior of W\ is similar to um) by using the su(2) (8> su(2)-coherent state. As 
a result, the obtained result is better in comparison with the exact first excitation energy. 

Further, we have improved the function p{K) by introducing the function p(K) in 
Eq. flTO]) . Here, we attach the suffix 2 to the quantities derived by using the function p(K) 
in Eqs. (jl8j) . In Fig. 3, the equilibrium energies are depicted in order to compare them with 
the exact ground state energy and E\ corresponding to the case p(K) = 1. The refine- 
ment by using the newly introduced function p(K) is also seen in the frequency uo in FigJH 
transparently. In the weak interaction strength region, the refinement is further obtained 
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in comparison with u\ with p(K) = 1. In the strong intearction strength region, the ob- 
tained frequency runs almost parallel to the exact one as a function of G. In addition to the 
above-mentioned refinement, even in the intermediate region of the interaction strength, the 
behavior is similar to the exact result and a fairly good result is obtained. 

In conclusion, the refined results for the first excitation energy around the equilibrium 
given by the coherent states have been obtained in the formulation of the idea in which the 
function g(K) reveals the similar asymptotic behavior to those of the su(2) ® su(l, 1)- and 
the su(2) (g> sw(2)-coherent states. The results shown in this paper teaches us that, following 
the increase of the interaction strength, the su(2) (g> su(l, l)-symmetry changes gradually to 
the su(2) (g) sw(2)-symmetry. An interesting future problem is to investigate which coherent 
state can describe the above-mentioned change of the symmetry. 
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